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Abstract 

We present an algorithm, based on the iteration of conformal maps, 
that produces independent samples of self-avoiding paths in the plane. 
It is a discrete process approximating radial Schramm-Loewner evolu- 
tion growing to infinity. We focus on the problem of reproducing the 
parametrization corresponding to that of lattice models, namely self- 
avoiding walks on the lattice, and we propose a strategy that gives rise 
to discrete paths where consecutive points lie an approximately constant 
distance apart from each other. This new method allows us to tackle two 
non-trivial features of self-avoiding walks that critically depend on the 
parametrization: the asphericity of a portion of chain and the correction- 
to-scaling exponent. 



1 Introduction 

The self-avoiding walk (SAW) is the prototypical lattice model for polymer 
behavior [Ij . It is defined as the uniform distribution over the set of all fixed- 
length nearest-neighbor walks on some lattice, such that no site is visited more 
than once. Self-avoidance reflects what are known as excluded volume effects in 
polymer science [2]. The universal aspects of the SAW have been the subject of 
study for decades in the physical and mathematical literature. While some of its 
properties are well understood and rigorously established (a thorough account 
can be found in (Hj), it still poses some difficult (and some seemingly impossible) 
problems. Mathematical and theoretical advances in the study of the SAW 
have always been paralleled by constant efforts for devising new algorithms and 
numerical strategies [4^^ 6 , and since it is one of the simplest non-trivial models 
it can serve as a test ground for novel algorithms in polymer science. 
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A connection has been studied in two dimensions — both numericaUy and 
analytically — between the critical SAW (in the limit where the number of steps 
goes to infinity) and a continuum model, called Schramm-Loewner evolution 
(SLE). SLE is a one-parameter family of stochastic processes in the complex 
plane producing random curves (traces) with conformal invariance "built in" . 
It has been conjectured that when the parameter is equal to | the scaling limit of 
self-avoiding walks is obtained |7] (this is the case we will be focusing on in this 
paper) . Later, numerical evidence has been given in favor of this correspondence, 
both in the half-plane [8] and in the whole-plane [9] geometries. 

What we propose here is an algorithm for sampling self-avoiding paths in 
the plane, based on a discretized version of SLE. Essentially, discrete paths are 
built by iterative composition of rotations together with one simple conformal 
map that takes a small circle and pulls a slit out of it. Thanks to the way 
SLE works, it is possible to efhciently produce independent samples, since the 
algorithm is based on simple Brownian motion, which is very easy to sample. 

Self-avoiding walks on the lattice have a natural parametrization, which 
corresponds to counting the number of steps along the walk. As long as one 
considers observables that depend only on the support of the walk, the corre- 
spondence with SLE curves is well-understood. But most of the quantities of 
interest in polymer physics do depend on the labeling of points along the chain 
and can not be matched with their SLE analogues, since SLE curves come with 
their own uncorrelated parametrization. Actually, the problem of finding a sen- 
sible definition of natural parametrization for SLE curves is still debated in the 
mathematical literature |10| lllj . From a numerical point of view, one needs an 
affordable way of generating SLE samples with the parametrization correspond- 
ing to the proper time of lattice models. One such method was introduced and 
studied by Kennedy [S] and will be briefly reviewed in Section [3] 

We hereby introduce a new method, based on the observation that the SAW 
— even when both the number of steps goes to infinity and the lattice spacing 
goes to zero — is such that the euclidean distance between two consecutive 
points on the chain is constant throughout the chain itself. We require the same 
property for the SLE discrete trace {7ti}, trying to attain an approximately 
constant step length 

|7n-7n-l|~A. (1) 

Discrete SLE chains are constructed by iterative composition of conformal maps, 
each one being responsible for producing a step 7n_i — >■ 7„. The method we 
propose for keeping an approximately constant step length does so by rescaling 
each step according to the Jacobian of the conformal map that acts on the 
corresponding segment. 

We focus on whole-plane SLE, so that our algorithm explicitly concerns 
the SAW in the plane, but the same reasoning and methods could be easily 
translated to other geometries, such as the chordal one. For instance, minimal 
modification is needed in order to treat ensembles of self- avoiding walks with 
fixed endpoints lying on the boundary of the domain. 

Section [2] is a brief introduction to the discrete process that we refer to 
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as discrete whole-plane SLE, which is the central object of interest lying at 
the heart of the algorithm; Section [3] introduces the issues about the choice of 
parametrization; Section |3] describes the numerical strategy used to reproduce 
the natural parametrization of lattice models in the framework of discrete SLE; 
numerical results for the SAW are presented in Section [SJ we measure the 
asphericity of an inner portion of SAW — which is a highly parametrization- 
dependent quantity — and we discuss the first correction-to-scaling exponent. 



2 Discrete whole-plane SLE 

We are not going to describe SLE here (the interested reader can find all the 
details in many excellent reviews, such as [121 HSl Ull HSl IS])- The aim of this 
section is to present a discrete process approximating radial SLE growing to 
infinity, i.e. a measure on curves with one end-point in z = 1 and the other at 
oo, living on the complex plane minus the unit disc, C \ D. A more in-depth 
presentation can be found in [S]. An analogous discrete process in the chordal 
geometry was introduced in [17] , where its convergence to SLE was also studied. 
Let us consider an ordered set of points 

7„eC\D, n^l,...,N. (2) 

We will call such a set trace or chain. We are going to describe a stochastic 
process whose outcomes are such traces. Let {6j} and {Aj} be two sequences 
of real numbers with Aj > 0. Consider the maps 

_ jz + 1)^ - 2c-^^z - (z + 1) ((z + 1)^ - 4c-^^zf 

These are conformal maps of D onto D\ [0|°'(1), l) , whose action can be described 
as growing a slit [0|^(1), l) inside D along the real axis, of length 



1 - 1 - e^^ (2 - e-^^ - 2Vl-e~^A . (4) 



cPj is the inverse map of the solution at time Aj to the Loewner equation in the 
disc 

f / \ . / Nexp(iat) -I- /t(z) 

dtft{z) = /t(z) -— — -, fo{z) = z, (5) 

exp(iat) - Jt[z) 

in the special case where the driving function is a constant at = 1. Notice that 
^ is well-defined also on the boundary of D. 

By complex inversion we can then define a family of maps growing slits in 
the complement of D in C: 
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^ - ^ cx 

Figure 1: (Color online) Composition of conformal maps giving rise to the 
trace. The first (top- left) arrow corresponds to the slit map the second one 
(clockwise) to the rotation _R„, the third one to 0„_i, the fourth one to Rn-i 
and so on. The last black arrow corresponds to the last rotation i?i . The dashed 
blue arrow is the complete map g^- Blue crosses in the last picture identify the 
points {7,1} that constitute the trace. 

The conformal maps (pj send C\D onto C\D minus a slit on the real line, whose 
length is controlled by Aj0. Let us construct a trace by intertwining such maps 
with rotations of the complex plane. We will consider the images of the point 
z = 1 under the chain of maps obtained by alternately composing rotations and 
slit mappings, as follows: 

^ Rl ° (f>l ° R2° (t>2° ■ ■ ■ ° Rn° </'n(l), (7) 

where 

Rj (z) = z exp (i(Sj) (8) 

are rotations whose angles are the parameters Sj. Let us call gn the n-th com- 
posed map, for use in the next sections: 

g„(z) = i?i O 01 O i?2 O 02 O • ■ ■ O -Rn ° (t)niz), (9) 

so that 7„ = gn(l)- In words, we traverse the sequences {Aj} and {Sj} back- 
wards from n to 1 and for each j we compose a slit mapping of parameter Aj 
with a rotation of angle Sj. Refer to Figure [TJ At first a single slit is grown, and 
the point z = 1 gets mapped onto the real axis, some distance away from the 
disc, namely at 0n(l)- Then the universe is rotated and another slit is grown 
by application of 4>n-i, so that the base of the previous slit — which still lies 
on the unit circle after the rotation — will be sent somewhere on the new slil|. 
Notice that in general the shape of the old slit gets distorted because of the 
action of (pn-i- The process goes on until the first map is reached; at that time, 
a chain of n points has been produced. Notice that, since the composition in 
^ goes backwards from n to 1, adding a step on the tip of the trace without 

^ See the discussion in Section |4] for the relation between Aj and the length of the slit. 

^ We suppose for clarity that the rotations are small enough for this to be the case, but even 
if the base of the slit remains on the circle the procedure explained here produces nonetheless 
similar sets of points {7n}. 
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changing the rest of the trace itself means inserting 0„+i at the rightmost place 
in d?]) and then recomputing the whole chairH. 

The trace obtained of course depends on the two sequences {Sj} and {Aj}, 
so that a measure on the latter induces a measure on the former. We will draw 
6j and in such a way that their relation be that of space versus time for the 
one-dimensional rescaled Brownian motion. The easiest way of doing so is to 
draw the Sj 's as Bernoulli variables in the set { y/ kAj , — y^wAj"} — where k is 
a positive real number — and we shall do so in the following. We still have the 
freedom to choose the "time steps" Aj. A careful choice of the latter is what 
allows us to reproduce the parametrization of the lattice models. 

The discrete process we have introduced here is expected to correspond to 
whole-plane SLE with parameter k in the appropriate limit (essentially, N — ^ 
oo), and in particular to the SAW when k = |. Notice that the presence of the 
unit disc as a forbidden region is expected to become irrelevant if one looks at 
the trace sufficiently far away from the origin, so it is not surprising that the 
lattice counterpart is full-plane SAW. 



3 The choice of parametrization 

The parameters {A^} defined in the previous section are the time steps of the 
discretization. They represent the time at which the Loewner evolution with 
constant driving function is to be evaluated in order to produce the slit maps 
defined in ([3]), which are the building blocks of the discrete evolution. If one 
is interested in reproducing the actual SLE, where time is indeed a continuous 
variable, one will want to have these parameters scale to zero. The choice of 
how they do so entails a choice of parametrization on the resulting object. For 
instance, taking a constant A^ = A and then sending A to zero yields SLE 
parametrized by capacity, which means that the curves have linearly increasing 
capacit}0. Different definitions can give rise to different parametrizations. This 
is not an issue when one is interested only in parametrization-independent fea- 
tures of the curves, such as the fractal dimension, the multifractal spectrum, 
the distance of the curve from some given point, or the probability of passing on 
the left or right of an obstacle, to name a few. But it becomes crucial when one 
focuses on parametrization-dependent observables of the lattice models. Such 
are for instance the distribution of a given point inside a SAW, the gyration 
tensor, properties related to the detailed shape of the walks and the universal 
quantities that describe the approach of an iV-step chain to = oo, such as 
the correction-to-scaling exponents (see Section [S|) . Our goal is to find a choice 
of {Aj } capable of reproducing the natural parametrization of the lattice walk 
models. 

In the following we will mostly use letter n to denote an index running from 1 to N 
(labeling the points on the chain) and letter j to denote an index running from n to 1 (labeling 
the incremental maps that build up the n-th composed map) , but this is not a strict distinction, 
since they actually label the same sequences. 

* The (logarithmic) capacity is defined as the logarithm of the coefficient of the term 
in the expansion around oo of gn{z). 
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The first thing one notices when producing discrete chains with the algorithm 
described in Section [2] is that the parametrization by capacity yields points that 
get further and further away from each other as the process goes on. The step 
size 

In = \ln - 7n-l| (10) 

diverges when n — >■ oo. One can then choose to scale the time steps in such 
a way as to compensate for this. It turns out (see [9] for more details) that 
choosing 

A,-r\ (11) 

at least definitively in j, cancels out the drift in the distance between consecutive 
points. Notice that this choice is non-random, meaning that the reparametri- 
zation scheme does not depend on the realization of the stochastic process. 
Instead, the values of {A^} are chosen before the actual simulation takes place. 
Unfortunately, this strategy does not give the correct parametrization (see also 
[T5] . where the same reasoning is applied to the half-plane case). For instance, 
as far as the spatial distribution of the fc-th point (for a given k) along the chain 
is concerned, it gives exactly the same results as the parametrization by capac- 
ity. What happens is that the scaling form ((TT]) only ensures that the average 
distance between consecutive points be constant (the average is over the 
realizations of the stochastic process). But fluctuations around this average still 
retain all their correlations, since we are reparametrizing in a naive, non-random 
fashion. 

One solution to this problem was proposed and studied by Kennedy [H [18] , 
and is particularly adapted to the case when one is interested in the position 
of just a single point along the curve. The idea is to grow the trace with its 
parametrization by capacity and stop the growth when a fixed "length" has been 
reached. A definition of length for a discretized fractal object can be introduced, 
which turns out to be naturally dependent on a fixed length A, that is the scale 
the fractal length (or variation) is measured at. Some care must be taken when 
using this method, because of the dependence of the results on the choice of the 
scale A. One would like to send A to zero, in order to measure the variation 
at a finer and finer mesh, but at the same time A can not become too small as 
compared to the step size of the discretized trace, otherwise wilder and wilder 
rounding problems would completely spoil the computation. Moreover, one 
wants to send the total length of the curve to infinity, which by scale invariance 
amounts to shrinking the unit disc down to a point, so as to approach the truly 
whole-plane geometry. One is therefore confronted with a tricky double limit, 
which increases the effort to be put into the analysis of the numerical data, and 
can blur the estimation of the errors. 

The method based on fractal variation is especially suited for producing a 
single point on the chain at a given value of the parametrization. The strategy 
we shall adopt here is different. We aim at producing a discrete trace where the 
step sizes /„ are strictly constant throughout the chain. Stopping the discrete 
growth after a fixed number of steps will then be automatically equivalent to 
choosing the stopping time when a fixed value of the fractal variation is reached. 
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The advantages of this strategy are manifold: one obtains an essentiahy arbi- 
trary number of points equally spaced in the natural parametrization at the 
same cost as producing only the last one, and no computationally-delicate dou- 
ble limit is present, so that no additional scaling analysis must be performed. 
Moreover, no prior knowledge of the fractal dimension d is needed. 



4 The numerical strategy 

A close relative to the slit mapping in dSj first appeared in the literature about 
diffusion-limited aggregation, or DLA. DLA — introduced by Witten and Sander 
|19j — is a kinetic model where finite-sized particles perform random walks (one 
at a time) from infinity until they stick irreversibly to a cluster, which grows 
from a seed placed at the origin. Hastings and Levitov [50] took advantage 
of the conformal symmetry inherent to this model and proposed an algorithm 
which turns out to be similar to what we use for simulating Schramm-Loewner 
evolutions. The algorithm works as follows. The seed of the growth is the unit 
disc. At each time-step, an angle 9j is chosen with the uniform distribution in 
[0, 27r). A conformal map 4>Aj,ej is applied, that creates a bump of fixed area 
Aj centered at e'^^ . Then, another 9 is chosen, the maps are composed, and so 
forth. Note that — as is the case for Equation (|9]) — if gn-i is the map that 
grows the cluster up to the (n— l)-th deposed grain, then the map that grows the 
cluster up to the n-th grain is obtained by first applying the incremental map 
4>A„,e„ and then <?«_!, which is to say that the incremental maps are composed 
in the opposite order than usual. 

This growth process satisfies an even stronger version of the domain Markov 
propert50, which is one of the crucial characteristics of SLE [21 , since now the 
growth of the cluster at a specific time does not even depend on where it last 
grew, so not only does the future not depend on the past, but — modulo a 
conformal transformation — it does not depend on the present either. 

An important technical aspect of this algorithm is that one wants to grow 
bumps of approximately equal size. But peripheral bumps have undergone sev- 
eral conformal maps and have thus changed their shape and size to a great 
amount. In general, by the time the whole cluster has been built, the n-th 
bump created (by 4>A,^,e„) has been subject to the action of gn-i- To compen- 
sate for this rescaling, one wants to create bumps of different sizes, depending 
on the whole history of maps they will be subject to in the remainder of the 
growth process. As a first approximation, as long as the new (n-th) bump is 
sufficiently small, it is natural to try and correct only for the Jacobian factor 
of the previous composed map, calculated at the place where the new 
bump is to be created, because this is the rescaling factor that will affect the 
shape of the bump at first order in its size. The n-th bump size should then be 



^ It is essentially the usual Markov property, but "up to" a (time- and realization- 
dependent) conformal transformation of the domain. 
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Since this strategy seems to give satisfactory results, it is very natural to try and 
apply it to the numerical reparametrization of SLE: rescaling the step sizes of 
the approximated SLE trace by the dilatation factor given by the Jacobian (very 
similar ideas were also fruitfully exploited in |22| , where multifractal spectra for 
Laplacian walks are computed). 

The size L„ of the n-th slit grown is a function of the time-like parameter 
A„ which controls the capacity of the incremental map at step n: 



A„ ^log \ , (13) 
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l + Lr. 

as can be seen by inverting (|4]) . One wants to rescale A„ so that L„ gets rescaled 
by a factor given by the Jacobian 

J„ = |.g;_i(e'*")|, (14) 

in analogy with (|12|) . 

Unfortunately, there happens to be a great obstruction to this program, due 
to the fact that SLE satisfies "only" domain Markov property, instead of the 
complete independence of DLA steps that we discussed above. If we look at 
Equation (fT2|) we see that rescaling the step-sizes is possible because of the 
independence of On (the space- like variable) from An (the time- like variable). 
This independence in DLA stems from the fact that the distribution of the 0„'s 
is flat on [0, 27r) and does not change, so that one can operatively choose every 
step 9o,0i, . . . before performing the composition of the corresponding maps. In 
SLE, on the contrary, despite the fact that the steps satisfy the domain Markov 
property, the increments Sn are drawn with a Bernoulli distribution from the set 
{VkA^, — v/kA^}, which does depend on time, since it depends explicitly on 
the time-like parameter A„; on the other hand the Jacobian needed to rescale 
A„ is to be evaluated at exp(i(5„). Therefore, the problem is that we do not 
really know where to compute the Jacobian, until we have actually computed 
it! This is ultimately related to the fact that SLE is driven by a non-trivial 
stochastic process, so that A„ and (5„ are intertwined. 

The reader can find a depiction of how the length of the slit L„ and the 
angle Sj are related in Figure [H which is a polar graph (for —tt < 5„ < tt) for 
the position of the tip of the slit as a function of the angle — as can be found 
by inverting (fT3|): 



2 - exp i-Sl/^) - 2v/l - exp i~5l/n) 

Thus, the main problem with the foregoing approach is that Sn and A„ de- 
pend on one another, so that one does not know where to compute the derivative. 
One way to overcome this problem is the following. Expand the derivative of 
Qn-i (the map that grows the hull at step n — 1) around its zero, which occurs 
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Figure 2: (Color online) The locus of the points Rn o 0n(l) (dashed blue lines), 
as Sn and A„ take on their allowed values. The red circle is the unit disc; the 
red arrows are examples of the possible slits grown. 

at z = 1, and evaluate it at the point exp(i(5„), which is the point where the 
n-th slit is going to be placed: 

5;_i(e'^") = 5;-i(l) + (e'^" - l) glUi^) + ■■■ (16) 

This expression is accurate when Sn is small. On the other hand, we also want 
to approximate the change in length of the slit by the value of the derivative at 
the base, silently assuming that it does not change much along the slit. This 
approximation is justified by the fact that, by (jlSp . Ln{Sn) is proportional to 
dn for Sn small. 

By expanding the exponential, taking the modulus, and remembering that 
g'n{l) = for every n — as can be seen by taking the derivative of ^ with 0„ 
given by ^ and ^ — one obtains the Jacobian (ITU) at order \Sn\ 

Jn-l - \Sn\WLl{l)\ ■ (17) 

We want to rescale the length L„ of the n-th slit by J„_i, so we rewrite the 
equation relating A„ and L„ (fT3|) by substituting 

\Sn\ |5"-l(l)| 

where A is the desired step length (as in ([IJ, which represents our goal), and by 
making use of the Brownian relation 

\Sn\ - (19) 
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We obtain an equation which (if solved) gives the time-step A„ producing both 
the correct rescahng, at first order in and the right relation with the space- 
step 5„: 

e--- = if2 , V^Vfl , ' , (20) 

The actual sign of (S„ = ±\/kA^ is to be chosen at random, according to the 
Bernoulli nature of 

Unfortunately, Equation (PH)) is transcendental, and can not be solved ex- 
plicitly. A little thinking shows that for |g"_i(l)| large one expects a small A„. 
In fact, (pO)) implies that the combination |ff^'_i(l)| v^kA^ be divergent when 
An — > 0. A crude approximation is then obtained by expanding the left hand 
side in powers of A„ around 0, the right hand side in |(7"_i(l)| -\/kA„ around 
00 and matching the two behaviors at first order. This is the best one can do, 
since higher orders would require solving algebraic equations of order greater 
than 4 in VA„. Among the solutions we choose the positive one: 

A„ = I (21) 

2\/k|5„-i(1)| 

Numerical solution to ^ shows that ^ is off by 20% when |5"_i(l)| = 10 
and by ~ 5% when |5"_i(l)| = 100, for A = 1 and k = |. The approximation 
works (see Section [S]) because the typical values of |5^i'_i(l)| involved are large. 

Of course, the foregoing method can be used only if one has an effective 
means of computing the main ingredient; |<?"_i(l)|. It turns out that there is 
such a way. Straightforward calculations show (the details are in the Appendix) 
that 

n-2 

l5"(i)l = l<(i)ink"-i-.(r.)|, (22) 

where Tj is defined as 

Tj = Rn-j O O Rn-j+l O (t>n-j + l O • • • O i?„ O 0„ (1) . (23) 

Equation is a closed formula for the Hessian, in terms of |(/)^(1)| — which 
only depends on A„, see ([32]) — and the function </>' — see (jSO]) and pT|) . 
Notice that the points rj=o,...,n-2 must already be computed by the routine 
that produces the n-th point on the trace 7„ = Ri o (pi o ■ ■ ■ o i?„ o (/>„(1) 
— as was explained in Section [5] — so that computing |5^'(1)| adds very little 
computational load. At the n-th step of the algorithm — i.e. when producing 
the n-th point along the chain — one can compute the factor |(7"(1)|, which will 
be needed for producing the (n-l-l)-th point, simply by multiplying the constant 
|(^^(1)| together with all factors |(/)^_i_j(rj)| obtained at each composition that 
is performed to compute 7„. This amounts to performing an operation taking 
a time 0(1) for each composition, which sums up to 0{n) for the n-th point 
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— that requires n compositions — and finally to 0{N'^) for a complete iV-step 
chain. 

Let us schematically sum up how our algorithm for building an TV-step chain 
works: 

1. Set the constant A (we shall always fix A = 1) 

2. Set n •(— 1 and Ai <- 1 

3. Compute i5„ as ±y/KA„ with a random sign 

4. Setup a temporary variable D •(— 10^(1)1 calculated as in (15^ 

5. Set j n and z 1 

6. Cycle on j for computing 7„: 

(a) Apply j-th incremental map and rotation z Rj o (j)j{z) 

(b) If j > 1 multiply D by calculated as in ^ and 

(c) If j > 1 decrease j by 1 and repeat step [S] 

7. Set 7„ ^ 2; 

8. Compute A„+i as in ((2T|) with |.g"(l)| given by D 

9. If n < increase n by 1 and go back to step [3] 

An example of a chain obtained by this method is presented in Figure |31 
where it is compared with a chain obtained by non-random rescaling as in 



5 Asphericity and corrections to scaling 

Given a chain {jn} — both for a walk on the lattice and for a discrete SLE 
trace — one can define its gyration tensor, which encodes useful information 
about the shape of the walk. Since the process we have defined grows a chain 
towards infinity, it can not be compared to a whole A^-step walk on the lattice, 
because the latter displays finite-chain corrections close to its tip. For this 
reason, we define the internal gyration tensor, by following the definition of 
gyration tensor that is used in polymer science, but by taking into account only 
the first M monomers: 

G"'W = ^E(7r-7;)(7f-7f), (24) 

where the superscripts a and /3 take values on the x and y coordinates of a 
lattice site or on the real and imaginary parts of a complex number. Let qi(M) 
and q2{M) be the two (real) eigenvalues of the (symmetric) gyration tensor. 



11 



Global rescaling 
Full algorithm 





Figure 3: (Color online) An example of a trace obtained by the full algorithm 
described in Section|4](in blue), compared to a trace obtained with simple global 
rescaling of the steps (in red). The sequences of signs in (t„ = ±-v/kA„ (that is, 
of left/right turns) are the same. 

These quantities are not universal, but some of their combinations are believed 
to be, such as the asphericity: 



which is a measure of how spherical the object is, being for perfectly spher- 
ical objects (for which the two eigenvalues are equal) and 1 when one of the 
eigenvalues is (as happens for objects lying on a line). 

The (critical) limit we are interested in is when the number of steps TV goes to 
infinity, because this is the limit where the SAW displays its universal behavior 
and where the cut-off at length 1 introduced by the disc-shaped forbidden region 
becomes irrelevant for the discrete SLE. On the SAW side, moreover, we will 
want to let Af — > oo, so as to avoid corrections to scaling, but in such a way 
as to have M/N — > 0, since we want to be looking at a portion deep inside the 
walk. 

We have simulated the self-avoiding walk on the square lattice using the 
pivot algorithm [23l [21] and we have measured the internal asphericity as a 
function of M/N for N = 100 000 and M = k ■ 1000 with /c = 1, 2, . . . , 99. The 
results are in Figure HJ the curve plotted is expected to be universal, and to our 
knowledge has never appeared in the literature. In order to obtain the value at 
M/N = Owe perform a fit of the form f{M/N) = y^gAW _^ a{M/NY on the 




(25) 
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Figure 4: (Color online) The internal asphericity of a SAW as a function of the 
fraction of walk considered. The error (not shown in figure) is « 10~* for each 
point (about the size of the red crosses), but this is to be taken cum grano salis, 
since different points on the plot are obtained from the same set of walks and 
are therefore not independent. 



first few values (namely M/N G (0,0.2]) in the grapf0, obtaining 

ylSAW ^ 0.51343(5) 
[6* = 0.92(5)]. 

The same measure (now with AI = N) is performed on ensembles of 20 000 
discrete SLE traces (k = f ) of lengths = 100, 200, 400, 1000, 2000, 5000, (and 
5000 traces of length 10 000) produced with the algorithm described in the 
previous sections. The results are in Table [Ttfl. The simulation for N = 400 took 
approximately one hour on a 2GHz Intel Duo processor. 

In general, for an A^-step chain one expects the following behavior for the 
expectation value of a global observable O: 



N 





^ + ... 




N 


A2 ^ 




bo 


+ 


h 








Co 




f- 









(26) 



where the leading behavior (given by the exponent po) is corrected by analytical 

^ The results show a strong stabihty as the upper cut-ofT in M/N is changed, up to around 
the middle of the chain, where they start to drift. 

^ We have also run simulations (for N = 100, 200, 400) using a numerical solution to I I20I1 , 
for the sake of testing the approximation in I I21II . The results agree with those in Table [T] 
[0.4956(17), 0.5030(17), 0.5069(17) for N = 100,200,400 respectively]. 
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Table 1: The asphericity of discrete SLE traces for k = |, computed for several 
lengths N. 

N y^SLE 

100 0.4957(17) 

200 0.5019(17) 

400 0.5071(17) 

1000 0.5097(18) 

2000 0.5115(17) 

5000 0.5122(18) 

10000 0.5124(37) 



(with integer exponents) and confluent corrections. The exponent^ (Ai < 
A2 < • • • ) are universal. The asphericity is expected to be a constant in the 
large- limit, so that — 0. The scaling form then suggests a fit of the 
form f{N) = yt^LE ^ bo/N^\ which yields 

A^^^ = 0.51351(54) 
Ai = 0.694(60). 

The asphericity is in perfect agreement with the one obtained for the SAW. 
Our result for Ai agrees with the theoretical value (Ai — ]^ ~ 0.6875) ob- 
tained by conformal-invariance methods for polymers in good solutions [25] . 
The evaluation of Ai for SAWs has been the subject of debate in the past 
decades. In fact, the rich structure of corrections in ([25)1 makes it difhcult to 
extract precise values from numerical data. There is now strong evidence for 
the absence of a leading term with exponent Ai = ji on the square lattice, 
the first non-null confluent contribution having exponent | [26| . Notice that 
this does not configure a violation of universality, since the amplitudes of the 
corrections are model-dependent. By including the first analytical correction 
(at next-to-leading order) and fixing Ai to its theoretical value, we find 

ylSLE ^ 0.51356(44) 
ai = 0.003 ± 0.362, 

where the amplitude ai of the term ai/N is compatible with 0. 



6 Conclusions and outlook 

We have studied a stochastic process in the complex plane, based on discrete 
Stochastic-Loewner evolution, which gives rise to chains {7n} with approxi- 
mately constant steps |7„ — 7n-i|. The purpose was to build an algorithm for 

* The symbols Ai, A2, . . . are those conventionally used for these quantities: they do not 
have anything to do with the Aj's used in previous sections. 
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exactly sampling self-avoiding paths, by correctly reproducing the parametriza- 
tion induced by the scaling limit of lattice models, namely self-avoiding walks. 
The method is based on iterative composition of conformal maps, where each 
map acts by building a radial slit out of the unit circle, which will eventually 
become one of the steps of the discrete path. Each step has to be rescaled ac- 
cording to the Jacobian of the map that evolves it. This program encounters 
some technical hindrances — essentially due to the fact that rescaling a step 
actually changes the Jacobian. We showed that an alternative approach is pos- 
sible, by keeping track of the second derivative of the map. It turns out that 
the Hessian can be effectively computed in the framework of iterated confor- 
mal maps, since it can be expressed as a simple function of quantities already 
computed by the algorithm. 

By exploiting the powerful correspondence existing between SLE and SAWs, 
the algorithm presented here produces completely independent samples of self- 
avoiding paths from the origin to infinity in the plane, whose parametrization 
is the desired one — that corresponding to SAWs — , and does so in an af- 
fordable way, with a complexity 0{N^) for A^-step chains. This allows us to 
study parametrization-dependent observables of the SAW such as the internal 
asphericity and the leading correction-to-scaling exponent, whose determination 
is considered a challenging problem in the numerical study of polymers. The 
results we obtain are very accurate. 

The analysis has been carried out in the whole-plane radial geometry, but 
very little should be changed in order to adapt it to the half-plane chordal case, 
or to other restricted geometries of interest in polymer science. 

Interesting questions remain open. It is still not clear whether there exists 
a way of reproducing finite-chain effects by the foregoing techniques. It would 
be useful for instance to produce the correct distribution of the end-point of a 
SAW, which is fixed to infinity in the present study. On the other hand, a great 
advance would be to translate this method to other classes of critical polymers, 
such as the 6 point where the collapsing transition poses even more difficult 
problems to Monte Carlo methods due to the attractive interactions and the 
consequently entangled shapes. 
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A Computing the second derivative of g^^ 

In this Appendix we give a formula for the modulus of the second derivative of 
gn, which is needed in a crucial step of the algorithm. 
The ?i-th composed map ^ can be written as 

gn{z) ^ gn-l O Rn O (t>n{z) (27) 
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in terms of the {n— l)-th map for n > 1. We recall that i?„ and 4>n are the 

n-th rotation and slit mapping respectively and that Rniz) implicitly depends 
on the parameter 5n, while (/)„ depends on A„. The second derivative of (|?7]) 
reads 

9n = K-l ° Rn O Cf>n) {K ° i^f 
+ (R'n ° <Pn) {(f>n) {g'n^l ° Rn O • 

The latter expression simplifies when evaluating its modulus at z = 1, because 
0^(1) = — a hallmark of the singularity of the Loewner map in 1 — and 
\R'n ° 'Pniz)] — 1- One obtains 



(28) 



lff;'(l)l = IC(l)l|5:.-l°i?n°0n(l)|. 

Computing the first factor |'?!>'^(1)| is just a matter of differentiating 
times with (fyj given by ([3]), which yields 



z2 mi/z)Y 



with 



(z - 1) e-^" + 1) y/{z + iy -Ae-^'^z - z-l) +2z 



2zV(-z + l)^-4e-'^"z 



and 



IC(1)I = ^e^" 



(29) 

two 

(30) 

(31) 
(32) 



VI -e-^ 

For the second factor \g!^_i o o (f)n{l)\, instead, a closed recursion can be 
found by noting that, by differentiating (p7|) . one has 



(33) 



(34) 



Iffn-ll = \9'n-2 ° Rn-1 ° (t>n-l\ | i O 1 1 

= \9'n-2 ° Rn-1 ° (t>n-l \ \<P'n-l \ : 

SO that 

Ig'n-l O RnO 0n(l)| = \(f)'n-l ° Rn O 0n(l)| ' 

|5ri-2 ° ° 4>n~l O Rn O (/)„(1)| . 

The seed of the recursion is given by the first slit grown 

.gi(z) = i?io0i(z), (35) 
so that finally from (1^ . ([55]) and the derivative of ([55]) one obtains 

o i?„ O 0„(1)| • 

\(l)'„_2 ° Rn-1 ° 4>n-l ° Rn ° (l>nii)\ ■ (36) 



|(/>1 O i?2 O (/.2 O i?3 O (/>3 O • • • O O O i?„ O 4>n{l)\ 
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or, in a more compact form, 

n-2 

l5;'(i)l = IC(i)ink"-i-.(r,)| (37) 

with Tj defined as in 
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